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Abstract 

We explain the application of a recently developed analytic continuation method to extract 
the electromagnetic transition form factors for the nucleon resonances (iV*) within a dynamical 
coupled-channel model of meson-baryon reactions. Illustrative results of the obtained iV* — > jN 
transition form factors, defined at the resonance pole positions on the complex energy plane, for 
the well isolated P33 and D13, and the complicated P\\ resonances are presented. A formula has 
been developed to give an unified representation of the effects due to the first two P\\ poles, which 
are near the irA threshold, but are on different Riemann sheets. We also find that a simple formula, 
with its parameters determined in the Laurent expansions of irN — > irN and 7./V — > irN amplitudes, 
can reproduce to a very large extent the exact solutions of the considered model at energies near the 
real parts of the extracted resonance positions. We indicate the differences between our results and 
those extracted from the approaches using the Breit-Wigner parametrization of resonant amplitudes 
to fit the data. 

PACS numbers: 13.75.Gx, 13.60.Le, 14.20. Gk 
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I. INTRODUCTION 



The spectrum and form factors of excited nucleons are fundamental quantities for in- 
vestigating the hadron structure within Quantum Chromodynamics (QCD). The excited 
nucleons are unstable and couple strongly to meson-baryon continuum states to form nu- 
cleon resonances (called collectively as N*) in irN and 7iV reactions. It is well known that 
resonances locate on the unphysical sheets of complex energy plane and thus their properties 
can only be extracted from the empirical partial-wave amplitudes by analytic continuation. 
Recently we have applied an analytic continuation method developed in Ref.JH to extract 
N* pole positions^ from nN elastic scattering amplitudes determined in a fit (3| (JLMS) 
within a dynamical coupled channel model|I| (EBAC-DCC) of meson-baryon reactions. 

The scattering amplitudes obtained from a dynamical coupled- channels model of meson- 
baryon reactions, such as the EBAC-DCC model as well as the models developed in Refs. 
0] , are not available in an analytic form. They are obtained numerically by solving coupled- 
channels integral equations with meson-exchange driving terms. Thus, the predicted am- 
plitudes can only be analytically continued to complex energy plane numerically with a 
careful account of the analytic structure of the considered scattering equations. Obviously, 
the method depends on the dynamical content of each model. For EBAC-DCC model, this 
has been developed in Ref.JH and established using several exactly soluble models. In this 
paper, we explain how this method is used to extract 7*iV — > N* transition form factors 
from the multipole amplitudes determined from extending the JLMS analysis to investigate 
7iV -> 7riV[l0| and N(e, e'7r)N^ reactions. 

The electromagnetic 7*iV — > N* transition form factors give information on the cur- 



rent and charge distributions of N* and N. It can be shown [121. 1 131 ] that a resonance 
state > with a complex energy Mr can be defined as an 'eigenstate' of Hamiltonian 
H\i/j§* >= Mr\^§* > with the outgoing boundary condition for its asymptotic wave func- 
tions. Therefore the 7*iV — > N* transition form factor is defined by the current matrix 
element < ip§*\J em \N > which can be extracted from the residue R^n^n of electromag- 
netic pion production amplitudes at the resonance poles. To extract R n N,-y*N, we need to 
evaluate the on-shell matrix elements of 7*iV — > nN amplitudes on the complex Riemann 
energy sheet. As will be discussed later, the analytic structure of the considered coupled- 
channels equations for getting these on-shell matrix elements is rather complex and must be 
dealt with carefully. In particular, we need to develop a formula to give an unified represen- 
tation of the first two Pu resonances which are near the 7rA threshold, but are on different 
Riemann sheets. 

To illustrate our approach, it is sufficient to only present results for the well isolated 
resonances in P33 and -D13 and the complex P\\ partial waves. With only three complex 
parameters determined in the Laurent expansion of each partial- wave amplitude at resonance 
pole position, we present a simple formula which can reproduce to a very large extent the 
exact solutions of the considered model at energies near the real parts of the extracted 
resonance positions. This finding agrees with what was reported in an analysis [Tij of 71 N 
scattering amplitude within the Jiilich model jsf. Here we show that this formula is also a 
good approximation for 'yN — > ttN amplitudes. Despite that this formula is similar to that 
used in the analysis [l5- 12 1 using the Breit-Wigner parametrization of resonant amplitudes 
to fit the data, we find no simple relation between two approaches. 

In section II, we will briefly review the analytic continuation method developed in Ref.[l| 
and explain how it is applied to evaluate the on-shell amplitudes of ttN, 7*iV — > nN transi- 
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tions. Section III is devoted to explaining how the determined residues are used to extract 
the elasticity r] e i of N* — > ttN decay and the 7*iV — > N* transition form factors at resonance 
poles. The results for Pn, P33, and P13 nucleon resonances are presented in section IV. A 
summary is given in section V. 

II. ANALYTIC CONTINUATION METHOD 

Within the formulation |ij for EBAC-DCC model, the partial wave amplitudes of meson- 
baryon reactions can be written as 

T^ a (p',p;E) = tf3 ja (p , ,p;E) + t* a (p , ,p;E), (1) 

where a, (3 represent the meson-baryon (MB) states 7A^, vrA^, rjN, pN, aN, 7rA, and 

tl a (p', P] E) = ^eAP f -iE)[G N *(E)} hJ T a ,(p;E) (2) 

with 

[G N l] hJ (E) = (E - m N: )5 t , - E y (tf) . (3) 

Here i,j denote the bare A^* states defined in the Hamiltonian. m^* are their masses. 
The first term (called meson-exchange amplitude from now on) in Eq.([T]) is defined by the 
following coupled-channels equation 

tp,a{p',P]E) = v j3> a(p',p)+ / dqq 2 ^2 v ^i(P^ ( liE)G 1 (q,E)t lja (q,p;E) (4) 

where vp jCt is defined by meson-exchange mechanisms, and Gy(q,E) is the propagator for 
channel 7. The dressed vertexes and the energy shifts of the second term in Eqs.Q-Q are 
defined by 

r ad (p;E) = r ad (p)+ / dgg 2 ^t ai7 (p',g;P)G 7 (g,P)r 7ij (g) (5) 

E(E) itj = f rfgg 2 ^r 7 ,(g)G 7 (g,P)f 7J (g), (6) 
j c 7 

where T a ^{j)) defines the coupling of the i-th bare A^* state to channel a. 

Because of the quadratic relation between energy E and momentum p, there are two 
energy sheets for each two-body channel: the physical (unphysical) sheet is identified with 
Im(p) > (<) for the stable two-particle channels. Thus the scattering amplitudes of 
an n-channels model are defined on a Riemann energy sheet which consists of 2 n sheets. 
For the EBAC-DCC model, defined by Eqs.([T])-([nD, each sheet can be defined by symbol 
(ZnN, z rjN, z-k-kNi -ZttAj z pNi z ctn) , where z a could be p or u representing the physical or un- 
physical sheets of channel a. Note that an acceptable reaction model can only have bound 
state poles and unitarity cuts on the physical sheet (pppppp). The sheets from all other 
possible combinations of u and p are called unphysical sheets on which the scattering am- 
plitude can have poles. We are however only interested in poles which have large effects on 
scattering observables and therefore they must be on the sheets which are near the {pppppp) 
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physical sheet. These poles are called resonance poles, and other poles are called shadow 



poles. It is known [18J, |19| that a shadow pole near the threshold of a channel can also 
have large effects on scattering observables and must also be considered in the search. As 
analyzed in Ref. [l| using several exactly soluble models, these poles are in most cases on 
sheets where the open(above threshold) meson-baryon channels are on unphysical u sheets 
and the closed(below threshold) channels are on physical sheet. In below we first recall how 
the analytic continuation method we had developed in Ref.[l| is used to search for such 
resonance poles within EBAC-DCC model. We then describe how it is used to extract the 
residues of the extracted resonance poles from on-shell amplitudes. 

Since v a> p and the bare vertex T a ^ are energy independent within the EBAC-DCC model 
, the analytic structure of the scattering amplitude defined above as a function of E is mainly 
determined by the Green functions G 7 (g,£'). Thus the key for selecting the amplitude on 
physical sheet or unphysical sheet is to take an appropriate path of momentum integration 
C in Eqs.([T])-(IHD according to the locations of the singularities of the meson-baryon Green 
functions G a (p, E) as E moves to complex plane. This can be done independently for each 
meson-baryon channel. For channel with stable particles such as vrA^ and r)N, the meson- 
baryon Green function is 

which has a pole at the on-shell momentum p defined by 

E = ^m 2 M +p 2 + ^/m 2 B +pl (8) 

As an example, let us consider the analytic continuation of the amplitude to the unphysical 
sheet of the MB channel when the energy E is above the threshold Re(E) > vtlb + rriM 
and Im(E) < 0. The on-shell momentum p for such a E is on the second and the fourth 
quadrant of the complex momentum plane. As Im(E) becomes more negative as illustrated 
in Fig. [H the on-shell momentum (open circle) moves into the fourth quadrant. The 
amplitude on the unphysical sheet can be obtained by deforming the path C into C\ so 
that the on-shell momentum does not cross the integration contour. For energy below the 
threshold for the MB channel (E < rris + rn^j), the on-shell momentum p su t> is on the axis 
of positive imaginary. As the energy moves into the region of Re(E) < rris + ttlm and 
Im(E) < 0, p su b moves to the second quadrant of complex p-plane and does not cross path 
Ci, as indicated by the dotted curves in FigJT] Hence the amplitudes on the physical sheet 
of MB channel for energy below MB threshold can also be obtained by taking the path C\. 

For the channels with unstable particle such as the ttA, as an example, the Green function 
is of the following form 



where 



£a(p,£) 



{TA,nN(q)} 2 q 2 dq 



c 3 E- E n {p) - [{E K {q) + E N (q)) 2 +p 2 ]V2- 

(10) 
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FIG. 1: The shift of the on-shell momentum (open circle/solid circle) of the two-particle Green 
function Eq. ([7]) as energy E moves from a real value above/below the threshold energy to a 
complex value with negative imaginary part. C{ is the integration path for calculating Eqs.(HD -([()]) 
amplitude for E on the unphysical Riemann sheet. 

The 7rA Green function Eq. (Q has a singularity at momentum p = p x , which satisfies 

E - E n {p x ) - E A (p x ) - Z A (p x , E) = 0. (11) 

Physically, this singularity corresponds to the irA two-body 'scattering state'. There is also 
a discontinuity of the nA Green function associated with the irirN cut in H A , as shown in 
the dashed line in Fig. El where p is denned by 

E = E„(p ) + [(m w + m N f +p%} 1/2 . (12) 

Therefore, for Re(E) > m B + m M , 2m 7T + m N , the integration contour C must be chosen to 
be below the irirN cut (dashed line) and the singularity p x , such as the contour C2 shown 
in FigfSJ for calculating amplitudes on the unphysical sheet. 

The singularity go of the integrand of Eq. ffl"Q]) depends on the spectator momentum p 

E - E n (p) = [(EM + E N (q )) 2 +P 2 } 1 / 2 . (13) 

Thus go moves along the dashed curve, illustrated in Figj3l when the momentum p varies 
along the path C2 of Figj2j To analytically continue Sa(p, E) to the unphysical sheet, the 
contour C3 of Eq. ( ITU]) must be below q . A possible contour C3 is the solid curve in Figj3j 
We emphasize here that we can deform the contour C only in the region where the 
potential v a> p(p' ,p) and the bare iV* vertex T MB N *(p) are analytic. The contours described 
above are chosen only from considering the singularities of MB and ttttN Green functions. 
Thus they must be further modified according to the analytic structure of the considered 
v a,pip' iP) an d Fmb,n*(p) to obtain the scattering amplitude in the momentum region of 
interest. This consideration is specially necessary when we need to get the on-shell amplitude 
for extracting the residues of the identified resonance poles. The residue of the amplitude at 
resonance pole is evaluated from the 'on-shell' matrix element, where the on-shell momenta 
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FIG. 2: Contour C2 for calculating Eqs.(|4|)-([6]) for E on the unphysical Riemann sheet with the 
unstable particle propagators, such as Eq.([9]) for irA channel. See the text for the explanations of 
the dashed line and the singularity p x . 
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FIG. 3: Contour C3 for calculating the A self energy Eq, (llOp on the unphysical Riemann sheet. 
Dashed curve is the singularity qo of the propagator in Eq. (1130 . which depends on the spectator 
momentum p on the contour C2 of Figj2j 



are defined as M R = E*{p™ N ) + E N (p° n N ) for 71 N channel and M R = q^ n + E N (q° n N ) 

with Q 2 = (q°%) 2 — (q^'n 1 ) 2 f° r 1*N channel. Since on-shell momentum are in general 
closer to the real axis than momentum on contour C, the analytic properties of the meson- 
exchange potential has to be examined. For example, the t-channel meson exchange potential 
v m'b> mb(p iP) °f tli e EBAC-DCC model has singularities at 

A 2 -(p-tf) 2 = (14) 

with A = Em'(p') — Em(p) or Eb>(p') — Eb{p)- The form of Tmb,n*(p) is chosen such that 
its singularity is at the pure imaginary momentum. Thus the contours have to be chosen to 
also avoid these singularities. As an example we show in FigJUthe singularities associated 
with the 7rA channel at E — (1357, — 76i) MeV. The dotted line for irnN cut and the circle 
shows px are the singularities from the Green's function, as discussed above. The most 
relevant singularity of the meson-exchange potential in our investigation of electromagnetic 
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FIG. 4: The contour (solid curve) for calculating electromagnetic matrix element, po and p x are 
the singularities shown in Figj2j The dashed-dot curve is the singularity of the pion-exchange 
-yN —7- 7rA matrix element at E = (1357, — 76z) MeV. 
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FIG. 5: The real (left) and the imaginary(right) parts of the half-off-shell matrix-elements of the 
non-resonant potential v^n^-kA for P\\ partial wave at E = (1357, — 76i), as functions of the real 
part (Re(p)) of off-shell momentum. The solid(dashed) curve is from the calculation along the 
path C 2 {C' 2 ) shown in FigJU 



pion production amplitude is due to the t-channel pion exchange of — > 7rA, which is 
shown as the dashed-dot curve. Thus the integration contour has to be modified to the 
solid curve Ci in FigJU This can be understood from Fig. [5] in which we see that the 
matrix element (dashed curves) of non-resonant potential iva^jv encounters the cut around 
Re(p) ~ 170MeV with the path C 2 , but it varies smoothly (solid curves) along the path C 2 . 



III. EXTRACTION OF TRANSITION FORM FACTORS 

To indicate the essential features of our approach more clearly, it is useful to first briefly 
describe how the resonance parameters are defined in the previous investigations. The 
scattering amplitude Fp jCl between any two channels a and /3 is related to the S'-matrix 
element by Fp a = (1 — Sp j0c )/2i. Within the rather general theoretical framework discussed 
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by, for example, Dalitz and Moorhouse|13l]. Taylor [20|, and McVoy[21|], Fp >a at energies 
near a resonance pole position Mr is parametrized as a sum of a pole term and a constant 
non-resonant contribution 

F p>a (E Mr) ~ + B Aa , (15) 

where i?g iQ , is the residue at the pole position Mr, and the non-resonant amplitude Bp iOL 
is an energy independent complex number. By the unitarity condition imposed on the 
full S matrix F a p(E) at E — > Mr, the non-resonant term B a ^ is written in terms of an 
non- resonant S-matrix S B , which is unitary by itself (S B S B ^ = 1) 

B ft « = (16) 

Then the pole term of Eq.f JT5|) is defined by the partial width T a and a phase <p a arising 
from the presence of the non-resonant term B a ^ 



Rp, a _ e^^/2e^^/r Q /2 



(17) 



E-Er E-Mr 
It is important to note that for the — > ttN amplitudes we are going to consider, 



e'^J^jv is the electromagnetic 'jN — > N* form factor which clearly must be a com- 
plex number when the the non- resonant term B^n^n is present and 7 at 7^ 0. We will 
see that our formula are consistent with these earlier investigations and will yield complex 
7iV — > N* form factors. Our main advance is to provide their interpretations in terms of 
dynamics defined within the EBAC-DCC model. 

Here we mention that by introducing appropriate energy-dependence of Ii7i(Mr), R a p, 
and B a> p, the expression Eq. (|15|) is used in practice to fit the experimental data. This is 
the origin of the commonly used Breit-Wigner parametrization of the amplitude in physical 



energy region. In some recent analysis |15Hl7j based on such a Breit-Wigner parametrization, 
the extracted 7*iV — > N* form factors are reported as real numbers. Clearly, this is rather 
different from what one can interpret from the above formula used in the earlier analysis fl3l. 



20, 21 



We now explain that within the EBAC-DCC model, it is straightforward to extract the 
resonance parameters Mr, B a ^ and R a ,/3 of Eq. (fT5|) by performing a Laurent expansion of 
the T-matrix defined Eqs.([T])-([6]). We need to find poles of scattering amplitudes T a ^. In 
principle the pole of the scattering amplitude can be found in the meson-exchange amplitude 
t and/or resonance amplitude t R of Eq. (TTJ). However as pointed out in Ref. [141] . a pole M x 
of the meson-exchange amplitude t does not survive as a pole of the full amplitude when we 
introduce coupling with bare N* states, since there is an exact cancellation between the pole 
contributions from t and t R at E = M x . Furthermore the non- resonant term at resonance 
pole t(E = Mr) is finite. Thus, the resonance poles of EBAC-DCC, or any model with bare 
N* states, can be found by only analyzing t R defined by Eq.(T2]). Consequently, we only need 
to explain how the residues of resonance poles are extracted from the term t R . 

The pole positions Mr of t R are found from the zeros of the determinant of TV* propagator 
defined by Eq.© 

A(E = Mr) = det[G N l (E = Mr)] = . (18) 
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The pole term of the N* Green function can be expressed as 

XiXj 



E-Mr' 



(19) 



where i,j denote the bare N* state in the free Hamiltonian and \% represents z-th 'bare' 
resonance component of the dressed iV* and satisfies 



(20) 



If there is only one bare N* state, with G N l(E) = 1/(E — m^* — £ {E)), it is easy to see 
that 



1 



X yfl- Z>(M R ) 

where H'(M R ) = [d'E/dE]E=M R - If we have two bare N* states, Eq. (l20p leads to 



(21) 



Xi 



X2 



M R -m N *-Z 22 (M R ) 



A' {M R ) 



Mr - m N * - S 22 (M i? 



-Xi 



(22) 
(23) 



where A' (Mr) = [dA/dE]E=M R can be evaluated using Eq. (fT8l) . 

Now it is straightforward to see how the residues Rp, a and non-resonant term Bp_ a of 
Eq. (fT5l) can be extracted from the amplitude Tp >a defined by Eq.(JTJ. First we note that at 
E near the resonance pole M R , the full amplitude defined by Eq.(JTJ can be written as 



Tf,, a (p?,p?;E^M R ) = tp, a (p?,p?-,M R )+tf >a (p?,p?-,E^M R 



(24) 



where p° n is the on-shell momentum of channel a; e.g. Mr = E n (p™ N ) + E^(p°^ N ) for the 
nN channel, and t^ fil (p on ,p on \ M R ) is finite, as explained above. By using Eq.([2]) for the 
definition of t^ a and Eq. (ll9j) for the pole term of N* propagator, we can perform Laurent 
expansion of the on-shell element of Eq. (|24H to obtain 



tUp7'PT;E^m r ) 



1 L a 



E-Mi 



R 



+ Bp tCt + Bl a (E-M R ) + .... 



where 



Y.x^.ApT.Mr). 



(25) 



(26) 



Here the dressed vertex T a j is defined by Eq.([S]). The terms Bp >a and B^ a in Eq. (l2"5"j) depend 
on the matrix elements of meson-exchange amplitude t of Eq.([T]) 



d 



B^ = t p , a (pf^ Mr) + j E [{E- MR)t« a ( P °f jP ° a n ; E) 



lR ( on m on 



E=M R 



(27) 



The term Bl a can be calculated, but is not relevant to our following discussions. 



Let us now consider Eg. (125!) for a = (3 = nN case. We need to relate the residue T n N^nN 
of its pole term to the residue of the nN elastic scattering amplitude F n N,nN defined by the 
standard notation 

F nN ,ME) - " [ M r -E 1e ^ Mr ' (28) 

where S n N,nN is the partial-wave S-matrix. In terms of the normalization of EBAC-DCC 
model S nNtnN (E) = 1 - 2i[np on E 7T (p on )E N (p on )/E}T nNtnN (p on ,p on ; E) , we find that (p on 
stands for p^ N ) 

F^nMMr) = _ n P 0nE "(P 0n J E ^ P °^ T nN<nN (P° n i V° n , M R ) (29) 

Mr 

Keeping only the pole term of Eq. (I25[) in evaluating the above equation and using the 
definition Eq. (1281) . we then obtain 

i(f> p on En (p on ) E n (p on ) - R R 

Re = K ^ 7riV nN ( 30 ) 

The nN elasticity of a resonance is then defined as 

* = (31) 

With the similar procedure, we can perform the Laurent expansion of 7*iV — > nN am- 
plitude to obtain 

V R (n on \V R (n on (1 2 \ 
T nN MP° n > <T; E->M R ) = KN[P J"t + B %N>lN + ■■■ (32) 



E — M, 



R 



where q on is the ^N on shell momentum defined by Mr = go + E]s[(q on ) and the momentum- 
transfer Q 2 = {q on ) 2 —qo- As discussed in section I, a nucleon resonance can be interpreted [l^. 
131 ] as an "eigenstate " of the Hamiltonian H\ip§, >= Mr\%1) r , >. Then from the spectral 



expansion of the Low Equation for reaction amplitude T(E) — H' + H' E ^ H H', where we 
have defined H' = H — H with H being the non-interacting free Hamiltonian, we have 

J-TTNaNVP ,<7 ,E-)-M R ) - — h •• (33) 

t, — Mr 

Obviously, we can see that < \H'\q on , Q 2 >=< ipfj*\J^{Q 2 )e^\N > is determined by the 
electromagnetic current operator J^(Q 2 ). It must be a complex number since the resonance 
wavefunction ip§, contains scattering states. Comparing Eqs. (l32l and (|33|) . we then interpret 
r R N (q on ,Q 2 ) as the N R — > 7*iV transition form factor. As seen in Eq. (fT9|) . the resonance 
consists of all bare N* components and hence we have 

< ^\J»(Q 2 )e,\N >=J2Xif jN ,N : (q° n \Q 2 ) (34) 



Using the normalizations defined in Ref . [ll| and following the definition originally intro- 
duced for the constituent quark model|22|, the usual 7*iV — > N* transition form factors are 
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TABLE I: The extracted resonance poles (ReMn, —ImMn) MeV and elasticity r] e (Eq. ([3"Tj) ) are 
compared with the values listed by PDG [23]. 





M R (EBAC-DCC) 


location 


Mr (PDG) 


7] e (EBAC-DCC) 


7] e (PDG) 


P33 


(1211, 50) 


(u-ppp-) 


(1209 - 1211 , 49 - 51) 


100% 


100 % 


D 13 


(1527, 58) 


(uuuupp) 


(1505 - 1515 , 52 - 60) 


65 % 


55 - 65 % 


Pn 


(1357, 76) 
(1364,106) 
(1820, 248) 


(upuupp) 
(upuppp) 
(uuuuup) 


(1350 - 1380, 80 - 110) 
(1670 - 1770, 40 - 190) 


49 % 
60 % 

8 % 


55 - 75 % 
10 - 20 % 



related to our extracted from factors by 

A 3/2 (Q 2 ) = CY,XjF% d (Q 2 , Mr, A 7 = 1,X N = -1/2) , 
3 

A 1/2 (Q 2 ) = CJ2Xj^ n /Q 2 ,M r ,X, = -1,X N = -1/2) 

3 

-- C'£xjr*. Ntj (Q t ,M R ,X, = 0,\ 1 



Sl/2(Q 2 



v/V 



-1/2) 



where Xn and A 7 are the helicities of the initial nucleon and photon, respectively, and 



c = E^_^ x j(2j + l)(27r)3(2g ) 



Air 



(35) 
(36) 
(37) 
(38) 

(39) 



where K = (M| - m? N )/(2M R ). 



IV. RESULTS AND DISCUSSIONS 

In this section, we illustrate our procedures by presenting the results for the pronounced 
resonances in P 33 , _D 13 and the complex P n partial waves. We also investigate the extent 
to which our results can be compared with those extracted from using Breit-Wigner form of 
resonant amplitudes to fit the data. 

Before we present our results for electromagnetic form factors, it is useful to first discuss 
our results from nN scattering amplitudes, which were briefly presented in Ref . [3, 0| • The 
extracted pole positions (Mr) and elasticities r\ e defined by Eq. (l31~|) for P33, D\ 3 and Pn are 
compared with the values from Particle Data Group [23| in Table [H We see that our results 
correspond well with PDG, while only one Pn near 1360 MeV is listed by Particle Data 
Group ( PDG) |23j ]. The extracted residues Re 1 ^, defined in Eq. (l30l . for nN amplitude are 
compared with some of the previous works in Table [III We see that the agreement in P33 
and P13 are excellent. However, we see that the residues of the Pn resonances extracted by 
four groups do not agree well while we agree well with GWU/VPI only for the resonance at 
1356 MeV. 

In Table [J we also indicate the location of each pole on Riemann energy sheet. Since 
we only search for poles in the region where the open (above threshold) channels are on 
unphysical u sheet and close channels (below threshold) on physical p sheets, as described in 
section II, the quantity deciding which sheet each resonance in Table [T]is on are the branching 
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TABLE II: The extracted irN residues Re" 1 ^ denned by Eq. (|30p are compared with several previous 
results. 





EBAC-DCC 




GWU-VPI[15j 




Cutkosky[25j 




Jiilich[14j 






R 


<t> 


R 


<P 


R 


<t> 


R 


<S> 


P 33 (1210) 


52 


-46 


52 


-47 


53 


-47 


47 


-37 


L>i 3 (1521) 


38 


7 


38 


-5 


35 


-12 


32 


-18 


Pn(1356) 


37 


-111 


38 


-98 


52 


-100 


48 


-64 


(1364) 


64 


-99 


86 


-46 










(1820) 


20 


-168 






9 


-167 







points for each channel, Within JLMS fit they are (1077, 1486, 1216, 1363 - 33z, 1703 - 
75z, 1906 - 323i) MeV for (irN, r]N, ttttN, it A, pN, aN ) , respectively. For example, the 
Pn pole at 1357 MeV ( 1364 MeV) is below (above) the ttA threshold 1363 MeV and is 
on upuuupp (upuppp) sheets since both poles are above irN and ttttN channels and below 
rjN, pN and aN channels. Thus their residues are very different although their positions are 
very close, since they are on different Riemann sheets. These two-poles structure near the 



7rA threshold are also found in the earlier analysis of VPIj24j and Cutkosky and Wang|25 
and the recent analysis by the GWU/VPl[l5| and Julich[l4j groups. 

Our results presented in Tables I and II suggest that the resonance parameters of the 
pronounced and well isolated resonance poles, such as P 33 (1210) and Di 3 (1527), are rather 
safely determined by the structure of the empirical partial wave amplitudes as far as the 
employed models have the correct analytic properties in the region not far from the physical 
region. On the other hand, the residues of poles near threshold are sensitive to the dynamical 
content of the models, as we have seen in the considered Pn case. 

We now turn to presenting our results for 7*iV — > N* form factors A\(Q 2 ) and S\(Q 2 ). 
We first observe that for the isolated resonances in P 33 and D\ 3 , Eq. (l25l) and Eq. (l2T|) for 
7*iV — > irN multipole amplitudes at E — > M R can be approximated as the following simple 
form 

T^N^NiE — > M R ) = B wN7 n — ^ 7rN, lj i (40) 

t, — M R 

where the complex constants are evaluated at resonance position E = Mr 

B^ N = t^ N {p on , q° n ; Mr) + M R )t^ N (p on , q on ; E)\ , (41) 



E=M R 



R wN ,y N = rK N (p™,M R )A x (Q 2 ,M R )/C, (42) 

where C is defined by Eq. (l39"]) . We observe that the expression Eq. pOl) . evaluated with all 
constants except E kept at their complex values at pole position M R , is a good approximation 
in the physical region of E near W R = Re(M R ). Similar good approximation is also for the 
ttN — > nN amplitudes, as also reported in Ref . [3] • Our findings are shown in FigsJ6] and [7] 
for the P 33 and P>i 3 partial waves, respectively. The determined constants B nN)J * N , R^n^n, 
B n N,TrN RnN,nN, and M R for each case in Figsj6] and [7] are presented in Table [LTTl 

We now note that the expression Eq. (f40l looks similar to the commonly used amplitude 
with a Breit-Wigner parametrization 



7riV ( ^ A (Q i E ) 

E — [W R - 
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TABLE III: Extracted resonance parameters. Rp a and Bp^ a are for the irN elastic scattering 
amplitude F^n^n and multipole amplitudes El±,Ml± of the pion photoproduction. 





M R (MeV) 




K N ^ N (W~ 3 fmMeV) B nNnN {l^fm) 


-F33 


1211 - 50i 


36.1 - 37.7i -0.43 + 0.13i 


M 1+ (3/2) -2728 + 1436i -7.43 - 3.86i 
#i+(3/2) 175 + 118i -3.49 + 1.511 




1527 - 58i 


37.6 + 4.9i 0.06 - 0.08i 


M 2 _(l/2p) -224 - 61.6i 1.01 - 0.44i 
E 2 -{l/2p) -437 - 368i 4.25 + 0.36i 




1100 1150 1200 1250 1300 1100 1150 1200 1250 1300 1100 1150 1200 1250 1300 

E(MeV) E(MeV) E(MeV) 

(a) (b) (c) 



FIG. 6: Energy dependence of the ttN amplitude(a) and the 77r Af 1+ (3/2)(b) and £i + (3/2)(c) 
amplitudes of P33 channel. The solid circle(triangle) shows real (imaginary) part of the amplitude 
calculated using Eq. (fTpjh The solid (dashed) curve shows real(imaginary) part of the amplitude 
of EBC-DCC model. 



where T tot (E = Wr) and T n N(E = Wr) are called the total width and partial decay width 
for irN channel, respectively, and A BW (Q 2 ,q,E) is assumed to be real numbers. The en- 
ergy dependence of these widths as well as the phase factor <J) BW (E) are parts of the as- 
sumptions in those analysis, which of course will influence how the non-resonant amplitude 
B^N^'NiQ 2 ] E) is adjusted to fit the data. 

Eqs.( l40l) and ( 1431) have similar structure, but they have important differences. First 
Eq. (1401) is evaluated at complex Mr and hence the on-shell momentum q on and p on are 
also complex. On the other hand, all energy and momentum variables in Eq. (l43l) are real 
numbers denned by the physical energy E. The non- resonant amplitude B^n^n m Eq. (l4TJj) 
is obtained from a coupled-channel calculation, while B^N^NiQ 2 ', E) in Eq. (l4Ul) is often 




FIG. 7: Energy dependence of the the 7riV amplitude(a) and the 77r M2_(l/2p)(b) and i?2-(l/2p)(c) 
amplitudes of D13 channel. The solid circle (triangle) shows real(imaginary) part of the amplitude 
calculated using Eq. (I15D . The solid (dashed) curve shows real(imaginary) part of the amplitude 
of EBC-DCC model. 
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calculated from tree-diagrams of phenomenological Lagrangian with unitarization using ttN 
amplitude. Thus it is difficult to see that the helicity amplitudes of 7*iV — > N* extracted 
from these two rather different approaches can be compared. 

The two-pole structure of Pu resonances near ttA threshold poses a problem in inter- 
preting our results for the 7*iV — > N* form factors A\(Q 2 ). We note that Eqs. (l40l) is valid 
for each of these two poles, but they are on different Riemann surfaces. Thus we need to 
find a parametrization which carries the sheet information in representing these two-pole 



contributions. Here we follow the approach of Refs. |26M28l| and a similar formula used in 
extracting meson resonances 19 , 29[ . 



We first use Eq. fl40l) to write the nN — > nN and 7iV — > nN scattering amplitudes on the 
n A physical(a = p) and unphysical(a=u) sheet as 

Tflipf^, E -+ M$) = -— %^ + flW (44 ) 

L — M R 

where a, /3 represent nN or ^N channels. All parameters R^ a , B^\ and Mr are obtained 
numerically from the amplitude as described in the previous section. The above two ampli- 
tudes with a = u,p can be combined by using the following unified representation 

TUpT'PT, E-,M R ) = _ ^ + ^^ + B + B i (45) 

tj - Mr - ■JPirA 



where p n A is the nA on-shell momentum p x determined by Eq.(ll). We require Tp >a 
Tp'^f 1 at p n a = P^/^f 1 ■ This requirement for a = nN, '-/N and /3 = nN determine; 



unknown complex numbers R, R l , Mr, B, B 1 and 7 from known parameters Rp a , B^ a and 
Mr\ Neglecting small contribution of R 1 and B 1 , we then obtain 

T^ a (pf,p™, E^Mr) = - f' a + B^ a (46) 

E- Mr- 7P-rrA 

where 

aAp) i\A u ) 

(47) 
(48) 
(49) 
(50) 

With = 49 - 68i MeV, M { r ] = (1359 - 76z) MeV and pJ2 = -65 + 86z MeV, M { ^ = 
(1357 - 76i) MeV, we have M R = (1364 - 105z)MeV, 7 = -0.146 + 0.062z and R^n^n = 
(-12-47z)MeV. The quantities R { ^% at Q 2 can be obtained from T% N T* N of Eq.(liD]) and 
hence R n N,jN can also be calculated from using Eqs.f l47|) -f l50|) . By interpreting R^n^n an d 
RnN,fN of Eq.(f46l) as the residues of a pole and using the procedures described above, we 
can then extract the electromagnetic helicity amplitudes A\(Q 2 ) and S\(Q 2 ). 



7 


Mr° ] ~ 

Jp) 
Pit A 


M ( r ] 

(u) 
PnA 




Mr 


= M$ - 


(P) 
IPnA 








- idp { :i/dE) 








(u) 
-PnA 


R/3,a 




- -ydpH/dE) 


PnA n l3,a- 



ldp ( :l/dE) 
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FIG. 8: Energy dependence of the PhttN scattering(a) and the Mi_(l/2p) 77r amplitude (b). The 
solid circle (triangle) shows real (imaginary) part of the amplitude calculated using Eq. (|46p . The 
solid (dashed) curve shows real (imaginary) part of the amplitude of EBC-DCC model. 



We have found that the unified formula Eq. (]46j) is a good approximation for both 7riV 
and 77r amplitudes if Eq.(44) is evaluated in the physical region where E is near Wr = 
Re(MR). This is shown in Fig. [8] for the considered Pn partial wave. Although Eq. fj46|) 
is close to the commonly used Breit-Wigner form of Eq. (j4*3j) . it is difficult to compare the 
extracted 7*iV — > N* helicity amplitude A\(Q 2 ) with those from previous analysis using 
Breit-Wigner parametrization, for the same reasons discussed above for the isolated P33 and 
-D13 resonances. 

We now present in Table HVl our results for the 7*iV — > N* for the P33, P13 and Pu reso- 
nances at Q 2 = photon point. As comparisons, we also list several previous results [30l-33l 



which were extracted from using the Briet-Wigner parametrization of resonant amplitude. 
Our results are complex numbers, as expected from expression Eqs.(35)-(37). Here we men- 
tion that a recent nucleon resonance analysis [34| also yields complex helicity amplitudes. 

We observe in Table [TV] that the real parts of our results for P 33 and D13 are in good 
agreement with the listed previous results. For the P33 case, this good agreement is perhaps 
related to the fact that the imaginary parts of our results for this pronounced resonance is 
much smaller than their real parts. For P 13 , a more detailed analysis is needed to under- 
stand this comparison since P 13 involves large ttN inelasticity and our results have large 
imaginary parts. For P n resonances, the real parts of our results (2c-bw) calculated from 
using the unified form Eq. (j4"6j) do not agree with the previous analysis using Breit-Wigner 
parametrization Eq. fj4*3|) . This is perhaps also related to the fact that our results for each 
pole near it A threshold have large imaginary parts, as also seen in Table IIVI 

For P33 we can use the standard relation 3] to evaluate the N-A magnetic transition 
form factor G* M in terms of helicity amplitudes. The real parts of our results (solid circles 
connected by solid curve) in FigJH] are in good agreement with the results (open circles with 
errors) from the previous analysis [3, 17] using Breit-Wigner parametrization. In the same 
figure, we also show that the imaginary parts (triangles connected by dotted line) of our 
results are much weaker. This observation further suggests that our results could be close 
to the results from analysis based on the Breit-Wigner parametrization only for the cases 
that the imaginary parts of our results are small. 

For the Pi3(1527) resonance, our results are shown in Fig. [TU1 As an example in see- 
ing the difficulty in comparing our results with those extracted from analysis using Breit- 
Wigner parametrization, we also show the the results (open circles with errors) from CLAS 
collaboration [l?}]. Qualitatively, CLAS analysis is based on the Eq. fl43p with the choice of 
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TABLE IV: The extracted jN — > N* helicity amplitudes (A\ i n 10 3 GeV 1 ' 2 ) are compared with 
previous results. 







"PR A C 
JitiAb 


Arnat|oUJ 


Aiirensioii 


Dugger[32] 


pianpieai 66 I 


P 33 (1210) 


Ay 2 


-265+19i 
-129+44i 


-243 ± 1 
-129 ± 1 


-256 ±3 
-137 ±5 


-258 ± 5 
-139 ±4 


-266.9 ±1.6 ±7.8 
-135.7 ± 1.3 ±3.7 


£>i 3 (1527) 


A3/2 
Ay 2 


171+9H 
-31+29i 


167 ± 5 
-20 ±7 


147 ± 10 

-38 ±3 


143 ± 2 
-28 ±2 




Pll(2cbw) 
(1356) 
(1364) 


A 1/2 

A/2 
Al/2 


-28+20i 
-13+20i 
-14+22i 


-63 ±5 




-51 ±2 





(f> BW = 0. Thus their Breit-Wigner amplitude become pure imaginary at E = Wr with 
Wr taken from PDG. As discussed in the beginning of this section, i.e. expression Eq.(17), 
the phase factor e l ^ E ^ is a necessary consequence of the presence of the non-resonance term 
-B7rjv,7*7v(p, Q 2 ; E) under the unitarity condition. This difference between the CLAS analysis 
and the previous analysis (l3l. 18 . 21 ] should be noted in interpreting their extracted 7*A^ — > 
form factors. 

Despite the differences between two different analysis, we observe that the real parts (solid 
circles connected by solid curve) of our A 3 / 2 and A1/2 shown in Fig. [10] are qualitatively 
similar to the CLAS data. The imaginary parts (solid triangle connected by dashed curve) 
of our results, which are smaller than the real parts but still appreciable, are also shown 
there. Since the longitudinal parts of the amplitudes could not be well determined with the 
available data, the large differences between our results and the CLAS data seen in Fig. [TU] 
are not very surprising. 

Our results for the three poles of P n listed in Table I are shown in Fig. [11] Similar to 
the results at photon point presented in Table III, their imaginary parts (solid triangles) 
are comparable or larger than the real parts(solid circles) in magnitudes. We note that the 



a 




0.5 



Q 2 [ GeV/c 2 ] 



1.5 



FIG. 9: The magnetic N-A (1232) transition form factor G* M (Q 2 ) defined in Ref.0j. F D = 
1/(1 + Q 2 /b 2 ) 2 with b 2 = 0.71 (GeV/c) 2 . The solid circles (solid triangles) are the real (imaginary) 
parts of our results. The other data points are from previous analysis [351] . 



1(3 




FIG. 10: The extracted 7A 



solid circles (solid triangles) are their the real (imaginary) parts 
collaboration[17 ] . 



JV(Di3(1527)) form factors (A l/2 (a),A 3/2 (h),S 1/2 {c)). The 

The data are from CLAS 




FIG. 11: The extracted 7A ->■ A*(1356)(a),A* (1364)(b) and 7V*(1820)(c) transition form factors. 
The solid circles (solid triangles) are their real (imaginary) parts. 



momentum dependence of the helicity amplitudes indicates that the structure of iV*(1356) 
and iV*(1364) is quite different from JV*(1820). 

It is perhaps more appropriate to interpret our results calculated from using the unified 
form Eq.(44) for the two poles near the ttA threshold as the values associated with the Roper 
A r *(1440) resonance listed by PDG. This results are shown in Fig. [121 We again see that its 
imaginary parts (dotted line) are comparable or larger than real parts (solid line) in most of 
the Q 2 region. Here we also see that the contribution ( dot-dashed lines) from the determined 
bare 7 AT — > N* strengths play an important role in changing the sign of the real part at Q 2 ~ 
0.4 (GeV/c) 2 . This sign change of the bare 7-/V — > N* form factor is seen in some relativistic 
constituent quark model calculations [36|, |37j. This suggests that our bare parameters can 
perhaps be interpreted in terms of hadron structure calculations excluding the meson-baryon 
coupled-channel effects which is determined by unitarity condition. Here we mention that 
our real parts are qualitatively similar to the results from CLAS collaboration. But it is not 
clear how to make connection between two results since we have very appreciable imaginary 
parts. 
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FIG. 12: Helicity amplitude Axi% v of P\\ resonance at ttA threshold extracted from using the 
unified representation Eq.(46). Solid (triangle) shows the real (imaginary) part of the helicity 
amplitude. Solid square (connected by dot-dashed line) shows contribution of the bare form factor. 



V. SUMMARY 

We have explained the application of a recently developed analytic continuation method 
to extract the electromagnetic transition form factors for the nucleon resonances (N*) within 
the EBAC-DCC model of meson-baryon reactions. We discuss in detail how the contours for 
solving the considered coupled-channels integral equations are chosen to find resonance poles 
Mr and their residues. The formula for determining the 7*iV — > N* transition form factors 
A\(Q 2 ) and S\(Q 2 ), defined on the complex Riemann energy sheet, from the extracted 
residues are presented. 

We have found that the resulting Laurent expansions of the ttN — >■ ttN and 7A — > irN 
amplitudes, evaluated in the physical energy region, can reproduce to a very large extent the 
exact solutions of EBAC-DCC model at energies near E = Re(Mn). A formula has been 
developed to give an unified representation of the effects due to the first two Pu resonances, 
which are near the ttA threshold, but are on different Riemann sheets. Illustrative results 
for the well isolated P33 and -D 13 , and the complicated Pu resonances are presented. 

We discuss the differences between our results and those extracted from the approaches 
using the Breit-Wigner parametrization of resonant amplitude to fit the data. We find that 
there is no simple connection between these two different approaches, despite that some of 
the real parts of our results and the results from Breit-Wigner analysis agree qualitatively 
when the imaginary parts of our results are much smaller. 

To conclude, we emphasize that our form factors are defined in a well-studied theoretical 



framework [131. Il8l. l21j within which a resonance is an "eigen state" of the Hamiltonian with 
the outgoing boundary condition for the asymptotic wavefunction of its decay channels. 
Thus the electromagnetic transition form factors defined by < ^^«|J em |A > , which can 
be extracted from the residues of resonance poles, must be complex, since the resonant 
wavefunction ipfj, contains scattering continuum. This must be accounted for in comparing 
our results with those from using the Breit-Wigner form to fit the data and any hadron 
structure calculations of N-N* transition form factors, such as those from relativistic quark 



models |36|, [37|, Dyson-Schwinger mo dels [3 81] . and LQCD[3 
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